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Abstract 

We present a calculation of the two-loop top-stop-gluino contributions to Higgs production via 
gluon fusion in the MSSM. By means of an asymptotic expansion in the heavy particle masses, we 
obtain explicit and compact analytic formulae that are valid when the Higgs and the top quark are 
lighter than stops and gluino, without assuming a specific hierarchy between the Higgs mass and 
the top mass. Being applicable to the heaviest Higgs scalar in a significant region of the MSSM 
parameter space, our results complement earlier ones obtained with a Taylor expansion in the 
Higgs mass, and can be easily implemented in computer codes to provide an efficient and accurate 
determination of the Higgs production cross section. 
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1 Introduction 



With the coming into operation of the Large Hadron Colhder (LHC), a new era has begun in the 
search for the Higgs boson(s). This search requires an accurate control of all the Higgs production 
and decay mechanisms, including the effects due to radiative corrections [IJ. At the LHC the main 
production mechanism for the Standard Model (SM) Higgs boson, Hsm, is the loop- induced gluon 
fusion mechanism [2], 5(7 — )• Hsm, where the coupling of the gluons to the Higgs is mediated by loops 
of colored fermions, primarily the top quark. The knowledge of this process in the SM includes the 
full next-to-leading order (NLO) QCD corrections O [U [5l [U El [8] ; the next-to-next-to-leading order 
(NNLO) QCD corrections [9] including finite top mass effects [10]; soft-gluon resummation effects 
pT] ; the first-order electroweak (EW) corrections [121 IISl E] ; estimates of the next-to-next-to-next- 
to-leading order (NNNLO) QCD corrections [15] and of the mixed QCD-EW corrections |16j . 

The Higgs sector of the Minimal Supersymmetric Standard Model (MSSM) consists of two SU (2) 
doublets. Hi and H2, whose relative contribution to electroweak symmetry breaking is determined by 
the ratio of vacuum expectation values of their neutral components, tan (3 = V2/V1. The spectrum of 
physical Higgs bosons is richer than in the SM, consisting of two neutral CP-even bosons, h and H, one 
neutral CP-odd boson, A, and two charged bosons, H^. The couplings of the MSSM Higgs bosons 
to matter fermions differ from those of the SM Higgs, and they can be considerably enhanced (or 
suppressed) depending on tan /3. As in the SM, gluon fusion is one of the most important production 
mechanisms for the neutral Higgs bosons, whose couplings to the gluons are mediated by top and 
bottom quarks and their supersymmetric partners, the stop and sbottom squarks. 

In the MSSM, the cross section for Higgs boson production in gluon fusion is currently known 
at the NLO. The contributions arising from diagrams with quarks and gluons, with full dependence 
on the Higgs and quark masses, can be obtained from the corresponding SM results [H O [6l [71 [8] 
with an appropriate rescaling of the Higgs-quark couplings. The contributions arising from diagrams 
with squarks and gluons were first computed under the approximation of vanishing Higgs mass in 
ref. [T7], and the full Higgs-mass dependence was included in later calculations [3 [71 [HI [18]. The 
contributions of two-loop diagrams involving top, stop and gluino to both scalar and pseudoscalar 
Higgs production were computed in the vanishing-Higgs-mass limit (VHML) in refs. [191 l2L)j . whose 
results were later confirmed and cast in a compact analytic form in refs. [2H I22j. Finally, first results 
for the NNLO contributions in the limit of vanishing Higgs mass and degenerate stop and gluino 
masses were presented in ref. [23]. 

The VHML can provide reasonably accurate results as long as the Higgs mass is well below the 
threshold for creation of the massive particles running in the loops. For the production of the lightest 
scalar Higgs, this condition does apply to the two-loop diagrams involving top, stop and gluino, 
but it obviously does not apply to the corresponding diagrams involving the bottom quark, whose 
contribution can be relevant for large values of tan /3. In turn, the masses of the heaviest scalar H and 
of the pseudoscalar A might very well approach (or exceed) the threshold for creation of top quarks 
or even of squarks. Unfortunately, retaining the full dependence on the Higgs mass in the quark- 
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squark-gluino contributions has proved a rather daunting task. A calculation based on a combination 
of analytic and numerical methods was presented in ref. [23] (see also ref. |25j). but neither explicit 
analytic results nor a public computer code have been made available so far. 

However, results from the first year of supersymmetry (SUSY) searches at the LHC (see, e.g., 
ref. |26) ) set preliminary lower bounds on the squark and gluino masses of the order of the TeV, 
albeit for specific models of SUSY breaking. This suggests that - if the MSSM is actually realized in 
nature - there might be wide regions of its parameter space in which all three of the neutral Higgs 
bosons are somewhat lighter than the squarks and the gluino. Approximate analytic results for the 
quark-squark-gluino contributions can be derived in this case. In particular, ref. [27] presented an 
approximate evaluation of the bottom-sbottom-gluino contributions to scalar production, based on 
an asymptotic expansion in the large supersymmetric masses that is valid up to and including terms 
of 0{ml/m'^), 0{mh/M) and 0{m'^/M'^), where rrifj, denotes a Higgs boson mass and M denotes a 
generic superparticle mass. An independent calculation of the quark-squark-gluino contributions to 
scalar production, restricted to the limit of zero squark mixing and degenerate superparticle masses, 
was also presented in ref. |28] . confirming the results of ref. [27] for the bottom contributions. More 
recently, ref. [22] presented an evaluation of the quark-squark-gluino contributions to pseudoscalar 
production that is also based on an asymptotic expansion in the large supersymmetric masses, but 
does not assume any hierarchy between the pseudoscalar mass and the quark mass, thus covering both 
the top-stop-gluino and bottom-sbottom-gluino cases. 

Exploiting the asymptotic-expansion techniques developed in refs. [27| and [22], we provide in this 
paper an evaluation of the two-loop top-stop-gluino contributions to Higgs-scalar production valid 
when the Higgs and the top quark are lighter than stops and gluino, without assuming a specific 
hierarchy between the Higgs mass and the top mass. In particular, we provide explicit and compact 
analytic formulae which include terms up to ©(m^/M^), 0{mf /M^) and ©(ml/M^). The results 
presented in this paper complement the earlier ones of ref. j21) . which, being obtained via a Taylor 
expansion in the Higgs mass, are not accurate for a Higgs mass comparable to (or greater than) the 
top mass, as might well be the case for the heaviest Higgs scalar of the MSSM. Our formulae can be 
easily implemented in computer code^, allowing for an efficient and accurate determination of the 
Higgs-boson production cross section in the MSSM. 

The paper is organized as follows: in section [2] we summarize general results on the form factors 
for Higgs boson production via gluon fusion in the MSSM. Section [3] contains our explicit results for 
the contributions arising from two-loop top-stop-gluino diagrams, as well as a discussion of the renor- 
malization conditions for the parameters in the top/stop sector. In section H] we compare numerically 
the results of our asymptotic expansion in the heavy masses with the results of a Taylor expansion 
in the Higgs mass, up to and including terms of O^rn^/m^) and ©(m^/M^), discussing the regions 
of applicability of the two different expansions and the effect of different renormalization conditions. 
Finally, in the last section we present our conclusions. 

^An implementation of the MSSM gluon- fusion cross section in the POWHEG framework was presented in ref. [29) . 
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2 Higgs boson production via gluon fusion in the MSSM 



In this section we recall for completeness some general results on Higgs boson production via gluon 
fusion in the MSSM. The leading-order (LO) partonic cross section for the gg ^ 4> process (with 
(p = h, H) reads 
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where is the muon decay constant, OsifJ^R) is the strong gauge coupling expressed in the MS 
renormalization scheme at the scale /x^, and Tifj) is the form factor for the coupling of the CP-even 
Higgs boson (p with two gluons, which we decompose in one- and two-loop parts as 



= nl' + fnf + o{ai) . (2) 

The form factors for the lightest and heaviest Higgs mass eigenstates can be decomposed as 
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T-Lh = Tp (—sin a + cos a ^2) 



Hh = Tp (cosa^i + sina?^2) 



(3) 



where Tp = 1/2 is a color factor, a is the mixing angle in the CP-even Higgs sector of the MSSM and 
T-Li (i = 1,2) are the form factors for the coupling of the neutral, CP-even component of the Higgs 
doublet Hi with two gluons. Focusing on the contributions involving the third-generation quarks and 
squarks, and exploiting the structure of the Higgs-quark-quark and Higgs-squark-squark couplings, we 
can write to all orders in the strong interactions [21] 
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In the equations above Af = 1/ sin/3 and A;, = 1/ cos /3. Also, /i is the higgsino mass parameter in the 
MSSM superpotential, Aq (for q = t,b) are the soft SUSY-breaking Higgs-squark-squark couplings and 
9q are the left-right squark mixing angles (here and thereafter we use the notation = simp, = 
cos(/3 for a generic angle (p). The functions Fg and Gg appearing in eqs. (HD and dSJ denote the 
contributions controlled by the third-generation Yukawa couplings, while Dg denotes the contribution 
controlled by the electroweak, D-term-induced Higgs-squark-squark couplings. The latter can be 
decomposed as 



Dg = fGg + C2e, 
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where I^g denotes the third component of the electroweak isospin of the quark q, Qg is the electric 
charge and 9w is the Weinberg angle. 

The form factors Hi can in turn be decomposed in one- and two- loop parts as in eq. ([2]). The 
one-loop parts, Ti}^, contain contributions from diagrams involving quarks (q) or squarks (qi). The 
functions entering Tij^ are 
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where = Amymf^, and the functions Qq and G\i2 read 
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The analytic continuations are obtained with the replacement m\ — )• m| -|- it . We remark that in 
the limit in which the Higgs boson mass is much smaller than the mass of the particle running in the 
loop, i.e. T » 1, the functions Qlf' and behave as 
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The two-loop parts of the form factors, Tif^, contain contributions from diagrams involving quarks, 
squarks, gluons and gluinos. We point the reader to, e.g., section 2 of ref. [21] for explicit formulae 
showing how Tif^ (or, equivalently, W^^) enter the total NLO cross section for Higgs boson production 
in hadronic collisions. In the next section we present our new evaluation of the top/stop contributions 
to T-Lf^, based on an asymptotic expansion in the stop and gluino masses. 



3 Two-loop contributions to the Higgs-production form factors 

In the case of the lightest Higgs boson h, the top/stop contributions to the two- loop form factor Tiff 
are well under control. Typically, the mass ratios between the lightest Higgs and the particles running 
in the loops allow for the evaluation of the relevant diagrams via a Taylor expansion in the Higgs mass, 
with the zero-order term in the series - for which ref. [21] provides explicit analytic formulae - already 
a good approximation to the full result. In the case of the heaviest Higgs boson H, on the other hand, 
the assumption that it is much lighter than the particles running in the loops is valid only in a limited 
portion of the MSSM parameter space. In particular, might very well sit around or above the 
threshold for the creation of a real top-quark pair in the loops, in which case - as found in ref. [22] for 
the pseudoscalar - a Taylor expansion in m% would certainly fail to approximate the correct result 
for the Higgs-production form factor. To address this possibility, we present in this section explicit 
analytic results for the two-loop top/stop contributions to the form factors Hf^ that include terms up 
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Figure 1: Examples of two- loop diagrams for gg ^ (p that do not involve gluinos. 
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Figure 2: Examples of two-loop diagrams for gg ^ (j) involving gluinos. 



to ©(m^/M^), 0{m\/M'^) and 0{m^/M'^), without assuming a specific hierarchy between the Higgs 
mass and the top mass. 

The top/stop contributions to Tif^ come from two- loop diagrams such as the ones depicted in figs.[T] 
and[2j In analogy to what was done in refs. |21[ I27j . we can decompose the functions F^^, Gf^, Ff^ 
and Gf^ entering the two- loop parts of eqs. ([4]) and ([5]) as 
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The various terms in eqs. (jl3p - (jl6p can be split in the contributions coming from diagrams with 
(s)top (s)quarks and gluons ((7, figs. [1^ and[l]3); with a quartic stop coupling (4t, fig. [1^); with top 
quarks, stop squarks and gluinos (^, figs. [2^ and [2b): 



[X = t,ti,t2,C2e^] 



(17) 



Furthermore, we remark that the term Yt entering eq. (|14p contains only contributions from diagrams 
with a Higgs-top coupling, figs.[T^ and[2^, therefore y/* = 0. On the other hand, the terms Y"^^ , Y^^ and 
Y^2 in eqs. p^ - (fT6]) contain only contributions from diagrams with a Higgs-stop coupling, figs. [TJ), 
[It, and [3). 
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3.1 Top-gluon, stop-gluon and four-stop contributions 

The top-gluon, stop-gluon and four-stop contributions to the terms in eq. ()17p can be extracted 
from the existing literature, and we collect them in this section for completeness. We assume that the 
parameters entering the one-loop parts of the form factors Tii in eqs. ^ and ([5]) are expressed in the 
DR renormalization scheme at the scale Q. 

The contribution to the term Yf arising from two- loop diagrams with top quarks and gluons (fig.[lJi) 
must be computed for arbitrary values of Tt = 4mj /rn^ . It reads: 



2m^Yf = Cf 
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where Cp = 4/3 and Ca = 3 are color factors, and exact expressions for T^j'.f^ 



, J- ana L/j^y2 



as functions of xt = — n — — rj -|- 1) are given in eqs. (2.12), (2.13) and (3.8) of ref. [7j, 

respectively. 

The contributions to the terms 1^^ , Yj^ ^'^d yg2^ arising from two-loop diagrams with stop squarks 
and gluons (fig. [T)o) and from diagrams with a quartic stop coupling (fig. [TJ;) can, to the accuracy 
required by our expansion, be computed in the limit of vanishing m^. They read [21] 
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The term is zero, while the terms Yf and Y^^ can be obtained by performing the substitutions 



-f-)- t2 in eqs. p9j) and (j20|) . respectively. 



3.2 Top-stop-gluino contributions 

In this section we present our original results for the asymptotic expansion of the top-stop-gluino 
contributions in the stop and gluino masses. We retain in our formulae only terms that contribute 
to the form factors T-Li up to 0{m'^ /M"^), 0{'m?^/M'^) or 0(m|/M^), where M denotes a generic 
superparticle mass. Again, we assume that the one-loop parts of Hi in eqs. and ([5]) are expressed 
in terms of DR-renormalized parameters evaluated at the scale Q. The top-stop-gluino contributions 
to the term Yt, arising from diagrams with a Higgs-top coupling (fig. [2^), read 
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where Xj = /m^ [i = 1,2), and 5mf^^^ denotes the SUSY contribution to the top self-energy, in 



units of as/i^-, expanded in powers of mt up to terms of 0{m^ 
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The terms IZi and 7^2 in eq. (j22]) coUect contributions suppressed by mt/M and mf/M^ , respectively: 
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We recall that the function Ql^2('^) defined in eq. pT]) . while 13{t) and K. {t) are defined as 
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Finally, -^{^2^'' can be expressed directly as a function of r in terms of the other three functions: 
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The top-stop-gluino contributions to the terms and Y^2 , arising from diagrams with a Higgs- 



stop coupUng (fig. [2)3), read 
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while Yf can be obtained by performing the substitutions xi — )• X2 and S20t — ?• —526* ^q. ([28|) . 

We remark that some care is required in order to properly include in the form factors Hi only 
terms up to 0{m'^/M'^), 0{m\/M'^) and 0{m^ /M"^). In particular, in the calculation of the function 
Ff^, eq. ()13p . we must use the full formulae for Yf and in eqs. ()28p and ()29p . respectively. On 

the other hand, in the calculation of the functions G^^, F^^ and G^^, eqs. p^ ~ p^ . we must retain 
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only the terms in the first two lines of eq. (j28p for Yf , and only the term in the first line of eq. (j29p 

for y| . 

As a first check of the correctness of our calculation, we verified that by taking the VHML (i.e., 



0) in the formulae presented in this section, which implies Qy^{Tt) — )• —4/3, /C (rt 
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we would obtain by expanding the VHML results of ref. [21] in powers of mt up to and including 
0{m1/M'^). It is also straightforward to check that, by performing the trivial replacement t — )• 6 in 
the formulae presented in this section and then dropping all terms that contribute to the form factors 
beyond ©(mf/m^), 0{mb/M) and ©(ml/M^), we recover the results of ref. [2T| for the bottom- 
sbottom-gluino contributions. To this effect it must be kept in mind that Q^^^^ti,) and K}^{Th) are of 
0(m2/m2), while B{n) = 2- \n{-ml/ml) + ©(mg/m^). 



3.3 On-shell renormalization scheme for the stop parameters 

If the parameters entering the one-loop part of the form factors are expressed in a renormalization 
scheme different from DR, the two-loop results presented in the previous section must be shifted in 
a way analogous to that described in section 3.2 of ref. [21j, to which we point the reader for details. 
Indeed, in our "on-shell" (OS) scheme we adopt the same prescriptions as in ref. |21] for the input 
parameters that are subject to 0[as) corrections: the top and stop masses are defined as the poles 
of the corresponding propagators; the counterterm of the stop mixing angle Ot is chosen as to cancel 
the anti-hermitian part of the stop wave-function renormalization matrix; the trilinear coupling At is 
treated as a derived quantity, related to the other parameters in the top/stop sector by 
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Some differences with respect to the treatment in ref. j2T] arise, however, due to the fact that 
the results presented in that paper were obtained in the VHML for arbitrary values of the top mass, 
while the results presented here are valid up to and including terms of ©(m^/M^) and ©(m^/Af^), 
without assuming a hierarchy between m,^ and rrif Defining the DR-OS shift for a generic parameter 
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to the order relevant to our calculation, the explicit expressions for the shifts in the stop masses and 
mixing angle read 
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where the analogous expression for Sm^ can be obtained by performing the substitutions ti -H- t2 and 



S29t ~^ ~^2et ill ei- (jST]) . We also define 5s2et = 2c2et ^^t and 5c2et 
mass reads 
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where the SUSY contribution 6mf^^^ was given in eq. 
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. Finally, the shift for the trilinear coupling 
At can be expressed in terms of the other shifts according to 



5A, 



2 



S29t 



2 
t2 



mt 



{At + 11 cot /3) . 



(34) 



If the one-loop form factors are evaluated in terms of OS parameters, the two-loop functions in 
eqs. (fT3ll - (fT6]l must be replaced by 
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(36) 

(37) 

(38) 



In addition, the two-loop form factor S^ts contributions originating from the shift in At: 

2m2 



mt S29t 
6S/3 



SAt 



1 



nir 



m? 

12 . 



+ 



15 



6 At 







V m5 

\ ti 


m^ / 

(2 / J 



(39) 



Once again, care is required in order to properly include in the form factors Tii only terms up to 
©(m^M^), ©(ml/M^) and ©(m^/M^). In particular: 

i) the shifts (5m? , (5m? and 56*4 must be computed up to 0{mt) in Ft^ , eq. ([35]) . while they must 
be truncated at order zero in mt in G^^, F^^ and G^^, eqs. ([5U|) - pS|) : 

zzj in F'I^ , eq. the first occurrence of 5mt must be computed up to 0{m1), while the second 
occurrence, in the term proportional to m?, must be truncated at order zero in m^; 
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in) in Gf^, eq. ()36p . the first occurrence of Srut must be truncated at 0{mt), while the second 
occurrence, in the term proportional to J-y^^^Tt), must be computed up to O(m^); 

iv) finally, in eq. (p9]) the first occurrence of 5 At must be computed up to 0{mt) by means of eq. 
while the second occurrence, in the term proportional to m^, must be truncated at 0{m~[^). 

Ref. [28j provides formulae for the two-loop SUSY contributions to the form factors for scalar 
production in gluon fusion, in the OS renormalization scheme, also based on an asymptotic expansion in 
the superparticle masses but restricted to the limit of zero squark mixing and degenerate superparticle 
masses. We checked that our OS results agree with those of ref. [28] in the simplified limit considered 
in that paper, after taking into account a difference in the overall normalization factor and the fact 
that ref. [28j employs the opposite convention for the sign of /i with respect to our eq. ([50]) . 

4 A numerical example 

We will now discuss a numerical evaluation of the two-loop SUSY contributions to the form factors 
for scalar Higgs production in a representative region of the MSSM parameter space. 

The SM parameters entering our calculation include the Z boson mass = 91.1876 GeV, the W 
boson mass ruw = 80.399 GeV and the strong coupling constant as{mz) = 0.118 [30]. For the pole 
mass of the top quark we take Mt = 173.2 GeV [3T]. For the relevant SUSY parameters we choose 

mg = m^/ = /i = 1 TeV , At = 2 TeV , m-g = 800 GeV , tan /3 = 5 , (40) 

where niQ and mu are the soft SUSY-breaking masses for the left and right stops, respectively. For a 
given value of the pseudoscalar mass m^, the scalar masses nih and niH and the mixing angle a are 
computed including the leading one-loop corrections of 0{at) and the leading two- loop corrections of 
Oiasat) M- 

In fig. Owe show the real part of the SUSY (i.e., all except top-gluon) contributions to the two-loop 
form factor for heaviest-Higgs production, H"^, as a function of rrijj. Since, as mentioned above, 
is not a free parameter in our calculation, its variation is obtained by varying between 100 GeV 
and 500 GeV. For simplicity, in the computation of the form factor we neglected the small D-term- 
induced electroweak contributions. The left plot in figure E] is obtained assuming that the parameters 
nit, Jn-tj; ^t2 ^* entering the one-loop part of the form factor, Tilf, are expressed in the DR 
renormalization scheme at the scale Q = 1 TeV. In this case we extract the DR top mass mt{Q) from 
the input value for the pole mass Mt by means of eq. (B2) of ref. [32], and we interpret the input 
parameters mg, mu and At in eq. (flO|) directly as running parameters evaluated at the scale Q. The 
right plot, on the other hand, is obtained assuming that the parameters mt, m^_^, m^^ and 9t entering 
are expressed in the OS scheme described in section 3.2 of ref. [21]. In this case we identify mt 
directly with the pole mass Mt, and we interpret the input parameters mq, mjj and At in eq. (j40p as 
the parameters that can be obtained by rotating the diagonal matrix of the physical stop masses by 
the "physical" angle 9t, defined through eq. (37) of ref. [21]. 
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DR scheme, SUSY only qS scheme, SUSY only 




nijj [GeV] mjj [GeV] 

Figure 3: Real part of the SUSY contributions to 'H'j^, plotted as a function of m^. The choice of 
SUSY parameters and the meaning of the different curves are explained in the text. The plot on the 
left refers to the DR scheme, while the plot on the right refers to the OS scheme. 



In each plot, the dashed (blue) line represents the result obtained in the VHML, as given in ref. [21j . 
while the solid (red) line represents the result computed at the first order of a Taylor expansion in m%, 
i.e. it includes the effect of terms of 0{m\/m1) and 0{m\/M'^) which were also computed in ref. [21] 
but proved too lengthy to be presented in analytic form. The dot-dashed (black) line represents instead 
the result of the asymptotic expansion in the superparticle masses derived in this paper. The latter 
is applicable when both nrit and rriu are smaller than the generic superparticle mass M, as is indeed 
the case here since M w 1 TeV, but it does not require any specific hierarchy between and m^. 

The comparison between the dashed and solid lines shows that, as rrijj increases, the effect of 
the terms of 0{m%/m^) and 0{m%/M'^) becomes more and more relevant, and the VHML does not 
provide an accurate approximation to Ti^. Furthermore, the comparison between the dot-dashed and 
solid lines shows that, even if the inclusion of the first-order terms pushes the validity of the Taylor 
expansion up to larger values of rrifj, the Taylor expansion fails anyway when ttih gets close to the 
threshold for the production of a real top-quark pair in the loops. In that case one can use the result 
of our asymptotic expansion in M, provided that the latter is still considerably larger than m^. 

A few additional comments are in order concerning the comparison between the left (DR) and right 
(OS) plots in fig. [3l There is no reason to expect the plots to look similar to each other, first of all 
because the difference between the values of H"^ in the two schemes is compensated for, up to higher- 
order terms, by a shift in the value of the one-loop form factor, Hjf, and also because the different 
interpretation of the input parameters in the two schemes means that, by using the numerical inputs in 
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eq. (j40p for both schemes, we are in fact considering two different points of the MSSM parameter space. 
This said, a striking difference between the two schemes is visible in the behavior of the asymptotic 
expansion (i.e., the dot-dashed hne) around the threshold for the production of a real top-quark pair 
in the loops. The fact that in the DR plot the threshold is located at a lower value of itih than in 
the OS plot is an artifact, due to lower value of the MSSM running top mass with respect to the 
pole top mass (indeed, for our choice of parameters mt{Q) = 144.3 GeV). The much sharper behavior 
around the threshold of the dot-dashed line in the DR plot, on the other hand, can be traced back 
to the contribution of the first term in the right-hand side of eq. (j22|) for . That term reflects 
the fact that the running top mass of the MSSM (i.e., including the stop-gluino contribution) is used 
in the top-quark contribution to "H^, and it is canceled out by the last term of eq. (|36p if the pole 
top mass (or, for that matter, the running top mass of the SM) is used instead. Indeed, we checked 
that, in a "mixed" renormalization scheme in which the stop contributions to T-L]^ are expressed in 
term of running parameters (including the MSSM running top mass) but the top-quark contribution 
is expressed in terms of the pole top mass, the qualitative behavior of the dot-dashed line around the 
threshold would be similar to the one in the OS plot. 

To conclude this discussion, we show in fig. [5] the real part of the SUSY contributions to the two- 
loop form factor for lightest-Higgs production, "H^^, as a function of the pseudoscalar mass m^, which 
is varied in the same range used to produce fig. [3l The meaning of the different curves is the same as 
in fig. [3l and again the left plot is obtained assuming that the parameters entering the one-loop form 
factor T-L]^ are expressed in the DR scheme, while the right plot is obtained assuming that they are 
expressed in the OS scheme. 

In the MSSM the mass of the lightest Higgs scalar h is bounded from above, and for large enough 
values of the pseudoscalar mass it becomes independent of m^, as do the couplings of h to the top 
quark and to the stops. Indeed, for the choice of SUSY parameters in eq. (jlUj) our crude 0{at + Oitag) 
calculation of the Higgs mass yields nih < 123.8 GeV in the DR plot and m/j < 122.5 GeV in the OS 
plot, and all the curves in fig.|3]become essentially fiat for > 250 GeV. Due to the relative smallness 
of nth no real-particle threshold is crossed, thus the result of the asymptotic expansion (dot-dashed 
line) is rather close to the result of the Taylor expansion at the first order in m| (solid line). 

However, a comparison between the left and right plots of fig. H] shows that in the DR calculation 
the VHML result (dashed line) provides a less-than-perfect approximation to Hf^-, while in the OS 
calculation the effect of the terms proportional to m| is small, and the VHML result essentially 
overlaps with the other two results. This difference between the two schemes can again be traced to 
the contribution of the first term in the right-hand side of eq. (j22p . i.e. to the choice of renormalization 
scheme for the top mass entering the top-quark contribution to T-L]!^ ■ Even in this case we checked 
that, in a "mixed" scheme in which the top-quark contribution to is expressed in terms of the 
pole top mass while the stop contributions are expressed in terms of running parameters, the VHML 
would provide as good an approximation to 'Hfj^ as it does in the full OS scheme. 
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Figure 4: Real part of the SUSY contributions to Hff, plotted as a function of m^. The choice of 
SUSY parameters and the meaning of the different curves are explained in the text. The plot on the 
left refers to the DR scheme, while the plot on the right refers to the OS scheme. 



5 Conclusions 

The calculation of the production cross section for the MSSM Higgs bosons is not quite as advanced 
as in the SM. Indeed, a full computation of the two-loop quark-squark-gluino contributions, valid 
for arbitrary values of all the relevant particle masses, has not been made publicly available so far. 
Moreover, the complexity of such a computation is going to be reflected in results that will probably 
be too lengthy and computer-time-consuming to be efficiently implemented in event generators. An 
alternative approach consists in deriving approximate analytic results that can be easily implemented 
in computer codes, and that are valid in specific regions of the MSSM parameter space such as, e.g., 
when the Higgs bosons are somewhat lighter than the squarks and the gluino. 

In this paper we presented a new calculation of the two-loop top-stop-gluino contributions to the 
form factors for Higgs scalar production in gluon fusion. We exploited techniques developed in our 
earlier computations of the production cross section for the MSSM Higgs bosons [271 ES] to obtain 
explicit and compact analytic results based on an asymptotic expansion in the heavy particle masses, 
up to and including terms of C'(m|/M^), 0{m1/M'^) and 0{m^/M'^). We compared our new results 
with the VHML results of ref. [21], as well as with the results of a Taylor expansion in the Higgs mass, 
up to and including terms of 0(m|/m^) and ©(m^/M^), and we discussed the regions of applicability 
of the different expansions. We also discussed the effect of choosing different renormalization schemes 
for the parameters in the top/stop sector. 
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From the example presented in section H] it appears that, in the case of the heaviest Higgs boson 
iJ, the use of our asymptotic expansion becomes mandatory when rrijj approaches (or crosses) the 
threshold for the production of a real top-quark pair in the loops. It also appears that choosing the 
OS scheme for the parameters in the top/stop sector leads to a milder behavior of the two-loop form 
factor T-L"^ around the threshold. In the case of the lightest Higgs boson /i, whose mass is bounded 
from above in the MSSM, we need not worry about thresholds. However, our discussion showed that 
in the DR scheme the VHML provides a worse approximation to than it does in the OS scheme. 

Finally, we remark that the results derived in this paper for the production cross section can be 
straightforwardly adapted to the NLO computation of the gluonic and photonic decay widths of the 
MSSM Higgs bosons, in analogy to what described in section 5 of ref. |21] . 
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